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Nonequilibrium Electron Interactions in Metal Films 
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Ultrafast relaxation dynamics of an athermal electron distribution is investigated in silver films 
using a femtosecond pump-probe technique with 18 fs pulses in off-resonant conditions. The results 
yield evidence for an increase with time of the electron-gas energy loss rate to the lattice and of the 
free electron damping during the early stages of the electron-gas thermalization. These effects are 
attributed to transient alterations of the electron average scattering processes due to the athermal 
nature of the electron gas, in agreement with numerical simulations. 

PACS numbers: 78.47.+p, 78.20.-e, 78.66.Bz 



The average properties of a free electron gas, such as 
the average electron scattering rate and the electron-gas 
energy losses, play a central role in the fundamental and 
technological properties of metallic and semiconductor 
systems. Although these are basically determined by the 
individual electron interactions, a collective behavior, re- 
lated to the electron temperature, is observed for a ther- 
malized electron gas. A drastically different behavior is 
expected for an athermal electron gas with thus strongly 
modified average electron properties. 

With the advance of femtosecond lasers, such transient 
athermal electron distribution can now be created and 
probed. In particular, in the case of noble metal systems, 
it has been shown that screening of the Coulomb interac- 
tion and the Pauli exclusion principle effect drastically re- 
duce the electron-electron scattering probability, leading 
to a slow internal thermalization of the electron gas on a 
few hundred femtosecond time scale [l|, 0, H, 0] ■ This lon- 
glasting nonequilibrium situation offers the unique pos- 
sibility of analyzing the interactions of an athermal free 
electron gas with its environment and their dynamic evo- 
lution during electron internal thermalization. 

Using silver as a model system, we have investigated 
the properties of a nonequilibrium electron gas using a 
femtosecond pump-probe technique. We show that off- 
resonant probing (i.e., far from the interband transitions) 
of the optical property changes of a metal film after ul- 
trafast perturbation of its electron distribution, permits 
selective investigation of the free electron-gas energy loss 
rate to the lattice and of the average electron scatter- 
ing rate. The results yield evidence for an increase of 
both of these rates as the initially athermal electron gas 
evolves to a thermalized distribution, due to the interplay 
between the different scattering processes for strongly 
nonequilibrium electron distributions. 

Measurements were performed in silver using a fem- 
tosecond pump-probe technique with near-infrared pulses 
to satisfy the off-resonant conditions (the threshold for 
interband transitions in silver is t&ln ~ 4 eV @). The 
electron distribution of an optically thin film is perturbed 
by intraband (free electron) absorption of the transform 
limited 18 fs infrared pulses delivered by a Ti: sapphire 



oscillator (hoj ~ 1.45 eV). Electrons are thus excited 
from below to above the Fermi energy, creating a strongly 
athermal distribution (Fig. [TJ inset) that internally ther- 
malizes by electron-electron scattering on a 500 fs time 
scale [H, S 0- The induced reflectivity, AR/R, and 
transmissivity, AT/T, changes are measured at the same 
wavelength using a time-delayed probe pulse. This off- 
resonant probing permits measurements of weak features 
associated with changes of the electron-gas average prop- 
erties. 

A conventional pump-probe setup has been used with 
differential and lock-in detection of the probe beam trans- 
mission and reflection changes. The two cross-polarized 
beams are focused over a focal spot of 30 mm in diam- 
eter on a 20 nm thick polycrystalline silver film grown 
on a sapphire substrate. The high repetition rate (80 
MHz) and stability of our femtosecond system permits 
noise levels for AT/T and AR/R measurements in the 
range. 
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FIG. 1: Transient transmission (solid line) and reflection (dashed 
line) changes measured in a 20 nm thick Ag film for infrared ex- 
citation and probing (fku ~ 1.45 eV) at To = 295 K. The pump 
fluence is 190 /xj/cm 2 . The inset shows the equilibrium electron 
distribution (/o, dotted line) and the perturbed one (/, solid line); 
Ep indicates the Fermi energy. 

The measured transient reflectivity, AR/R, and trans- 
missivity, AT/T, changes are shown in Fig.[T]for a pump 
fluence of 190 /xJ/cm 2 that would correspond to a max- 
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FIG. 2: Short-time behavior of the normalized AT/T (solid 
line) and AR/R (dashed line). Also shown are the pump-probe 
cross correlation (dash-dotted line) and its integral (dotted line). 

imum electron temperature rise AT e ~ 90 K for a 
thermalized electron gas. Similar results were obtained 
for pump fluences in the range 500 — 4 /iJ/cm 2 (i.e., 
AT e ~ 200 — 2 K) . The observed changes show very dif- 
ferent short time delay behaviors. The AT/T rise closely 
follows energy injection in the electron gas [i.e., the in- 
tegral of the pump- probe correlation (Fig. [3J)] . In con- 
trast, AR/R exhibits a delayed rise and reaches its max- 
imum value after about 120 fs, indicating that AR/R 
and AT/T are sensitive to different electronic property 
changes. No modification of the measured responses was 
observed with silver film aging, indicating that oxidation 
does not significantly influence our results. 

In the weak perturbation regime investigated here, 
AR/R and AT/T are linear combinations of the changes 
of the real, Aei, and imaginary, Ae2, parts of the metal 
dielectric function. The coefficients of these combinations 
can be calculated from the known optical properties of 
silver Q taking into account Fabry-Perot effects [||. For 
our probe wavelength, AT/T essentially reflects Aei (the 
contribution of Aei to AT/T is about 5 times larger 
than that of Ae2) while AR/R reflects both Aei and 
Ae2, whose transient behaviors can thus be determined 
(Fig. [3]). A striking difference here is that, as AT/T, Aei 
exhibits no rise time while, as AR/R, As2 reaches its 
maximum value after about 120 fs. This demonstrates 
the different physical origins of the induced changes of 



the dielectric constant dispersive and absorptive parts. 

The metal dielectric constant e = £\ + ie-i can be writ- 
ten as the sum of free (intraband) and bound (interband) 
electron contributions: 

e{uj) = e b {w) - . Up . . , (1) 

where oj p is the plasmon frequency and 7 is the free elec- 
tron optical scattering rate. In noble metals the inter- 
band term e b is dominated by transitions from the full d 
bands to the Fermi surface states. Energy injection into 
the electron gas results in a spreading of the electron dis- 
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FIG. 3: Induced changes of the real part (Aei, solid line) and 
imaginary part (Ae2, dotted line) of the metal dielectric constant 
obtained from the measured transmission and reflection changes in 
the 20-nm-thick film. 



tribution around the Fermi energy, Mlp, and thus to an 
induced (reduced) absorption for interband transitions 
to electron states below (above) Kip. These transitions, 
however, take place at large photon energies and hence, 
for sufficiently small pump and probe photon energy (i.e., 
2huj < hflib), AE2 is only due to changes of the intraband 
term (Ae^). In contrast, Aei, which is sensitive to the 
integrated absorption spectrum, is dominated by changes 
of the interband term Ae\ (the change of the real part of 
the free electron contribution Ae™'™ is negligible; see be- 
low). This can be calculated using the Kramers-Kronig 
relation: 



Ae b iH oc 



(cu' -n F + O lf) )- 1 g(^)[/ (^) - fjhio')] 

(«' -n F + n lb ) 2 ~uj 2 

1 



duo 1 , 



(2) 



where undispersed d bands have been assumed, /o and 
/ are the equilibrium and transient electron distribution 
functions and g is the conduction band density of states. 



In off-resonant conditions, one can show that Ae\ is pro- 
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FIG. 4: (a) Measured transient change of the real part Aei of 
the dielectric constant on a logarithmic scale. The dashed line 
corresponds to an exponential decay with r = 900 fs. (b) Calcu- 
lated evolution of the real part change of the interband contribution 
to the dielectric constant, Ae\ (solid line), and of the excess en- 
ergy stored in the electron gas, AE e (dotted line). 

portional to the electron-gas excess energy, AE e 

Ae^(cj) oc J u;' 3/2 [f(hiu') - f (Hu)')]du' oc AE e . (3) 

The short time delay increase of Aei thus essentially re- 
flects the time evolution of the excess energy stored in the 
electron gas and is expected to rise instantaneously, in 
agreement with the experimental results. Furthermore, 
its delay yields information on energy transfer to the lat- 
tice permitting direct investigation of its dependence on 
the electron distribution. 

The temporal dependence of Aei is shown in Fig. [4] 
on a logarithmic scale after subtraction of its final back- 
ground value. The initial decay of Aei clearly shows that 
the electron-gas energy losses to the lattice are initially 
slow and increase over a time scale of a few hundred fem- 
toseconds (Fig. 0|). For longer times, a constant rate is 
reached, corresponding to an exponential decay (Fig. HJ, 
in agreement with the theoretical results for a thermal- 
ized electron gas (two-temperature model [1, Q) in the 
weak perturbation regime [3j . The time constant of the 
long term exponential decay yields an effective electron- 
phonon coupling time of ~900 fs, consistent with pre- 
vious measurements A similar behavior has been 
observed in gold films with, however, an additional fast 
transient peak due to transient induced interband ab- 
sorption at the probe photon energy (in contrast to sil- 
ver, 2hio > HVlib in gold and Ae^w) ^ for very short 
time delays). 

This increase with time of the effective electron gas- 
lattice coupling is a direct consequence of the noninstan- 



taneous internal thermalization of the electron gas. Dur- 
ing the excitation process, a very small number of elec- 
trons gains a large excess energy as compared to fee To. 
As a crude approximation, separating the electron gas 
into unperturbed and nonequilibrium electrons, only the 
latter ones can lose energy by phonon emission. Their 
number increases with time as electron-electron scatter- 
ing redistributes energy among the carriers: this leads 
to an overall increase of the energy loss rate to the lat- 
tice during the early stages of the internal electron-gas 
thermalization 0, Q ■ As internal thermalization is ap- 
proached, the above separation is no longer valid and 
a constant energy loss rate is eventually reached, corre- 
sponding to collective electron-gas-lattice interaction as 
described by the two-temperature model. This evolution 
from a quasi-individual to a collective electron behavior 
is responsible for the observed short time nonexponential 
excess energy decay. 

The above description of the interplay between electron 
thermalization and energy losses can be made more quan- 
titative by computing the electron-gas relaxation dynam- 
ics using the Boltzmann equation for electrons, including 
electron-electron and electron-phonon scattering. Nu- 
merical simulations were performed for a free electron 
gas using statically screened electron-electron Coulomb 
interaction Deformation potential electron-phonon 
interaction has been introduced using the Debye model 
for the phonon dispersion. The only parameter is the 
amplitude of the deformation potential coupling which 
has been set by imposing the computed long time (ex- 
ponential) signal decay to be identical to the measured 
one. The time dependent changes of the interband term 
Aei are calculated from the electron distribution using 
the band structure model of Rosei et al. [j| . 

The results are in very good agreement with the exper- 
imental ones, showing in particular the same transient 
evolution of the signal decay to an exponential behavior 
(Fig. QJ. As expected, this behavior is very similar to 
that of the excess energy stored in the electron gas with 
only a slight short time delay deviation (Fig. [4j. This is 
a consequence of the initial perturbation of the occupa- 
tion number of electron states far below the Fermi energy 
that leads to more resonant transient contributions to e\ 
[Eq. (J5J)]. Similar calculations performed assuming in- 
stantaneous internal thermalization of the electron gas 
show only an exponential decay of Ae^, in agreement 
with the two-temperature model. 

The conduction band electron states probed from the 
d bands are not perturbed and, in contrast to Aei, Ae2 
reflects only the modification of the intraband optical 
absorption. Using the Drude model, this is connected to 
the alteration of the average optical scattering rate 7 of 
the electrons: 

Ae 2 ^ f u 2 -i 2 \ A7 ^ A7 

£2 \UJ 2 +"f 2 J 7 7 
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Our results thus demonstrate a noninstantaneous in- 
crease of 7 with energy injection in the electron gas 
(Fig. [3]) , and thus the modification of the efficiency of the 
scattering processes with the electron distribution char- 
acter (athermal or thermal). Ae 2 , and thus A7, subse- 
quently decays and reaches a plateau with a temporal 
behavior comparable to that measured for Aei. 

For short time delays all of the excess energy is stored 
in the electron gas, and the increase of 7 can thus be at- 
tributed to changes of the electron scattering process effi- 
ciency as the electron distribution broadens in the vicin- 
ity of the Fermi energy (weakening of the Pauli exclusion 
effects). Although at room temperature electron-phonon 
interaction is the dominant process, the other scattering 
mechanisms (i.e., electron-defect and umklapp electron- 
electron scattering) can also contribute to the observed 
changes. Their computation in the general situation of 
a noncquilibrium electron distribution is far beyond the 
scope of this paper, and we have only performed crude es- 
timates using the results for thermal systems [lo|, EH • For 
our probe photon energy, the dominant changes of the op- 
tical electron scattering rate is expected to originate from 
the temperature dependent part of the electron-electron 
scattering rate [lol ]. For an electron temperature rise 
of 90 K and a relative contribution of electron-electron 
scattering to the total scattering rate 7 e _ e /7 ~5x 1CT 2 
A7/7 is estimated to be of the order of 1CT 3 , cor- 
responding to Ae 2 ~ 4 x 1CT 4 [Eq. with e 2 ~ 0.5 
0| ]. This is consistent with the measured value when 
one takes into account the approximations in comput- 
ing umklapp electron-electron scattering and the possible 
role of the interfaces. The observed Ae 2 rise time sug- 
gests that, as for electron-lattice energy exchanges, the 
electron-electron collision rate also evolves during energy 
redistribution in the electron gas with its probability in- 
creasing on a 100 fs time scale. 

With the lattice heat capacity being much larger than 
the electron one, as electrons thermalize with the lattice 
their excess energy strongly decreases, eventually mak- 
ing A7 e _ e negligible. The concomitant rise of the lat- 
tice temperature, AT^, leads to a significant increase of 
the electron-phonon scattering probability which is pro- 
portional to the phonon occupation numbers and thus 
to Ti,. This effect eventually dominates the electron 
scattering rate change, 7 e - p h finally reaching a constant 
value determined by AT L : A7 e _ ph /7 e _ ph ~ AT L /T L . 
Neglecting long time scale heat diffusion, the electron- 
lattice equilibrium temperature rise is estimated to be 
AT e = ATl ~ 0.7 K, corresponding to an estimated long 
term e 2 change Ae 2 ~ 1.1 x 10~ 3 , in very good agreement 
with the experimental value (Fig. [3j . 

Modification of the free electron scattering might also 
alter the real part of the dielectric function [Eq. |T])] 
with a contribution proportional to Ae 2 : Ae I 1 ntra (w) w 
27/wAe 2 (o>). With the probe photon frequency being 



much larger than the electron scattering rate (7 ~ 20 
meV Q), A£ I 1 ntra estimated from the measured Ae 2 is 
thus much smaller than the observed Aei, indicating a 
negligible intraband contribution to Aei. 

In conclusion, transient optical reflectivity and trans- 
missivity changes have been investigated in silver films 
using a femtosecond pump-probe technique with 18 fs res- 
olution. Experiments were performed in the weak pertur- 
bation regime in off-resonant conditions, permitting one 
to connect the deduced changes of the real and imagi- 
nary parts of the metal dielectric constant to fundamen- 
tal properties of the free electron gas, i.e., respectively, 
the electron-gas energy losses and the electron optical 
scattering rate. 

The results demonstrate that the electron-gas energy 
losses to the lattice increase with time during the early 
stages of the electron-gas internal thermalization, in 
very good agreement with numerical simulations of the 
free electron dynamics. This behavior reflects evolution 
from an individual to a collective electron-lattice type 
of coupling, mediated by electron-electron interactions. 
Induced free electron optical absorption has also been 
demonstrated with a rise time of about 100 fs, show- 
ing a noninstantaneous increase of the average electron 
scattering rate. This indicates a strong dependence of 
the probability of the electron scattering on the electron 
distribution, and theoretical investigations of these pro- 
cesses for off-equilibrium distributions would be partic- 
ularly interesting here. The induced scattering is likely 
dominated by the increase of electron-electron scattering 
for short time delays and electron-phonon scattering due 
to the lattice temperature rise as the full system reaches 
thermal equilibrium. 
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